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ABSTRACT 


The  objective  of  this  project  was  to  study  the 
feasibility  of  drying  grain  at  ambient  air  conditions  but  at 
pressures  below  atmospheric  and  compare  this  to  the 
traditional  ambient  air  drying  system.  The  study  method  was 
to  calculate  the  energy  requirements  per  unit  of  moisture 
removed  from  the  grain  by  the  two  systems.  The  dry  and 
wet-bulb  temperatures  prevailing  near  Edmonton,  Alberta, 
were  used  for  the  study.  The  results  indicated  that  the 
energy  requirements  per  Kilogram  of  moisture  removed  are 
less  for  pressures  below  atmospheric  than  above,  and  that 
the  energy  requirements  for  pressures  above  atmospheric 
increase  at  a  rate  greater  than  the  pressures  below 
atmospheric  when  either  the  grain  depth  or  airflow  is 
increased.  The  number  of  hours  available  each  day  to  dry 
grain  to  14%  moisture  content  were  more  and  the  number  of 
days  required  to  dry  wheat  were  less  using  pressures  below 
atmospheric  for  the  usual  drying  period  in  Western  Canada. 
Further  research  should  continue  that  could  lead  to  higher 
energy  efficiency  in  grain  drying. 
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1.  INTRODUCTION 


In  his  recent  book,  Hail  (1980)  estimates  the  annual 
loss  of  grain  after  harvest  to  be  10%  of  the  potential 
yield,  and  the  loss  of  hay  to  be  28%.  These  estimates  are 
for  the  United  States;  for  many  other  countries  the  losses 
would  exceed  these  figures.  Moisture  control,  primarily  by 
drying,  provides  an  opportunity  to  reduce  post-harvest 
losses . 

Grain  drying  has  become  common  across  the  Canadian 
prairies  during  recent  years  (Friesen  1978).  Instead  of 
drying  only  during  wet  harvest  seasons,  many  farmers  now  use 
a  drier  as  part  of  their  normal  grain  harvesting  system. 
There  are  a  number  of  reasons  for  making  drying  an  integral 
part  of  the  harvesting  system.  Some  of  these  include  earlier 
harvesting  by  starting  at  higher  moisture  contents  and  thus 
reducing  field  losses  due  to  shattering,  wind  and  insects, 
and  longer  harvest  days  by  using  extra  hours  of  harvesting 
in  the  morning  and  evening.  It  also  permits  the  use  of 
late-maturing  cultivars  and  a  more  efficient  use  of  land  and 
labor.  The  grain  quality  produced  can  be  superior  which  is 
worth  more  to  both  the  farmer  and  the  consumer. 

The  initial  cost  of  the  grain  drier  is  one  of  the  major 
deterrents  to  grain  drying,  but  where  large  amounts  of  grain 
are  dried,  fuel  efficiency  is  a  more  important  factor.  With 
the  increasing  interest  in  drying,  such  as  cereal  grains, 
and  the  increasing  fuel  prices,  there  is  a  need  for 
investigations  aimed  at  energy  conservation  in  drying.  This 
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study  was  undertaken;  therefore,  to  see  if  there  is  a  drying 
system  that  requires  less  energy. 


2.  LITERATURE  REVIEW 


The  traditional  method  of  drying  farm  crops  is 
ventilating  the  grain  with  heated  or  unheated  air.  Based  on 
a  1973  estimate,  3.8  million  Kiloliters  (1005  million 
gallons)  of  LPG  (liquid  propane  gas)  was  used  for  drying 
crops  with  heated  air  (CAST  1975).  LPG  is  the  most  popular 
fossil  fuel  for  grain  drying  in  Canada  and  the  United 
States.  The  ever - i ncreas i ng  fuel  prices  and  possible 
shortages  in  future  has  made  it  necessary  to  use  alternative 
fuels,  to  improve  the  efficiency  of  existing  systems,  and  to 
develop  new  systems  for  energy  conservation.  The  following 
is  a  brief  review  of  the  drying  systems  including  the  drying 
theory  based  on  vapor  pressure  differential. 


2.1  DRYING  THEORY 

Drying  takes  place  when  the  vapor  pressure  within  a 
hygroscopic  material  such  as  grain  is  greater  than  the  vapor 
pressure  of  the  surrounding  medium  (which  is  usually  air). 
The  equation  presently  used  in  calculating  water  vapor 
transmission  through  materials  is  based  on  a  form  of  Fick's 
Law  and  is  as  follows  (ASHRAE  1977)  : 

w  =  -  u  (dp/dx )  ( 1 ) 

Where  w  =  weight  of  vapor  transmitted  through 
a  unit  area  in  unit  time, 
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p  =  vapor  pressure, 
x  =  distance  along  the  flow  path, 

(dp/dx)  =  vapor  pressure  gradient,  and 
u  =  permeabi 1 i ty . 

As  shown  by  equation  (1),  if  the  vapor  pressure  in  the 
grain  is  higher  than  the  vapor  pressure  in  the  surrounding 
air,  the  transfer  of  moisture  would  be  from  grain  to  air 
thus  effecting  drying.  Conversely,  if  the  vapor  pressure  in 
the  air  is  higher  than  the  vapor  pressure  in  the  grain  the 
transfer  of  moisture  would  be  from  the  air  to  the  grain  thus 
increasing  the  grain  moisture. 

There  are  two  methods  of  increasing  the  vapor  pressure 
differential  between  the  grain  and  the  air.  Traditionally  in 
grain  drying,  the  differential  is  increased  by  heating  the 
air.  As  the  hot  air  comes  in  contact  with  the  grain,  the 
latter  is  heated  thereby  increasing  its  vapor  pressure.  As 
will  be  shown  later,  the  energy  requirements  of  this  method 
of  drying  are  significant,  and  in  industrialized  nations, 
this  method  is  largely  dependent  on  fossil  fuels.  The  other 
method  of  increasing  the  vapor  pressure  differential  is  used 
in  food  drying.  Food  processors  use  reduced  air  pressures  to 
increase  the  vapor  pressure  differential.  According  to 
Dalton's  Law,  reducing  the  total  air  pressure  will  reduce 
the  partial  pressure  of  the  dry  air  and  the  water  vapor,  and 
with  a  reduction  in  the  latter,  a  vapor  pressure 
differential  will  be  increased  (Todd  1981).  This  process  has 
also  high  energy  requirements  but  low  pressure  drying  is 
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used  to  dry  food  because  certain  foods  retain  their 
characteristics  of  taste  and  nutritional  quality  better  when 
dried  under  low  pressure  and  low  temperature  than  at  high 
temperature.  However,  the  drying  apparatus  for  creating  low 
pressure  is  expensive.  Agricultural  crops  are  not  as 
sensitive  to  high  temperature  drying,  so  the  less  costly 
equipment  is  often  selected  by  many  farmers. 

Researchers  interested  in  grain  drying  have  attempted 
to  find  a  system  that  requires  less  energy  than  the  heated 
air  system  traditionally  used  by  the  farmers.  One  system  is 
based  on  ambient  air  where  the  air  is  pumped  through  the 
grain  mass.  Drying  will  occur  because  there  is  usually  a 
small  vapor  pressure  differential  between  the  grain  and  the 
ambient  air.  This  system  is  generally  considered  efficient 
but  it  has  some  other  limitations  which  are  discussed  later 
in  this  chapter . 

A  number  of  researchers  have  studied  the  energy 
requirements  of  a  variety  of  grain  drying  systems.  The 
results  of  six  most  commonly  used  drying  systems  are 
discussed  below  in  order  to  identify  the  typical  energy 
requirements  for  grain  drying. 


2.2  HIGH  TEMPERATURE  DRYING 

High  temperature  driers  use  a  large  flow  of  hot  air  to 
dry  the  grain  in  a  period  from  less  than  an  hour  to  several 
hours.  Air  temperatures  range  from  10°C  above  ambient  to 
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80°C,  the  upper  limit  being  determined  by  the  grain  and  its 
intended  use  (Hall  and  Davis  1979). 

The  energy  requirements  for  high  temperature  driers 
vary  with  drier  design,  drying  temperature  and  airflow.  A 
typical  energy  requirement  is  about  5000  KJ/kg  of  moisture 
removed  (Hall  1980).  The  statistic  includes  the  energy  to 
heat  the  air  and  drive  the  fan  with  the  fan  energy 
representing  about  10%  of  the  total. 

The  quality  of  grain  dried  with  heated  air  is  often 
lower  than  that  dried  without  heat  (Foster  1973).  The 
degradation  of  the  grain  quality  characteristics  such  as 
viability  or  milling  properties  is  a  function  of  the  air 
temperature.  Higher  air  temperatures  reduce  the  grain 
quality.  Because  of  the  reduced  grain  quality,  and  the  high 
energy  requirements  with  high  temperature  drying,  other 
systems  of  drying  are  frequently  used. 

2.3  LOW  TEMPERATURE  AND  AMBIENT  AIR  DRYING 

Low  temperature  drying  is  a  method  of  drying  grain  when 
small  amounts  of  heat  are  added  to  increase  the  temperature 
a  few  degrees  above  the  ambient.  When  no  heat  is  added,  it 
is  called  ambient  air  drying.  Research  has  shown  that  low 
temperature  and  ambient  air  drying  maintains  the  grain 
quality,  and  the  energy  requirements  relative  to 
conventional  dryinf  techniques  are  low  (Gustafson  et  al. 

1976;  Morey  et  al .  1976). 


i.. 


Otten  et  al.  (1977)  performed  some  low  temperature 
drying  experiments  near  Guelph,  Ontario,  using  corn.  They 
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reported  a  total  energy  input  for  heating  the  air  and 
driving  the  fan  from  3700  to  4400  kd/kg.  The  main  reason  for 
the  variation  in  the  energy  requirements  was  due  to  the 
higher  airflow  rates  and  a  variation  in  ambient  temperature 
(-17.5  to  9.5°C)  during  the  drying  period. 

Morey  et  al.  (1978A)  conducted  some  drying  experiments 
at  three  different  locations  (St.  Cloud,  Minnesota;  Des 
Moines,  Iowa;  and  Indianapolis,  Indiana)  to  evaluate  several 
fan  management  strategies  for  ambient  air  drying  of  corn. 

The  energy  requirements  (fan  only)  varied  from  900  to  2300 
kd/kg  of  water  removed.  In  a  similar  study  Morey  et  al. 
(1979)  performed  some  experiments  to  evaluate  the 
feasibility  of  ambient  drying  of  wheat.  They  found  the 
energy  requirements  (fan  only)  were  1700,  2000,  and  2400 
kd/kg  for  Fargo,  Grand  Forks,  and  St.  Cloud,  Minnesota, 
respectively.  These  energy  requirements  are  less  than  that 
reported  for  low  temperature  drying  and  much  less  than  the 
commonly  used  conventional  hot  air  driers.  Also  the  minimum 
energy  requirements  are  less  than  the  latent  heat  of 
vaporization  for  grain  moisture  (2800  kd/kg)  as  reported  by 
Brooker  et  al.  (1974).  This  is  because  the  heat  energy  for 
vaporization  comes  largely  from  the  sensible  heat  of  the 
ambient  air. 

Despite  the  success  obtained  with  low  temperature  and 
ambient  air  drying,  it  can  be  risky  (Shove  1973A).  The 
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wetter  the  grain,  the  greater  the  risk  of  spoilage.  Any 
system  for  drying  with  little  or  no  added  heat  should, 
therefore,  be  based  on  an  understanding  of  the  relationship 
between  the  temperature  and  the  grain  moisture  content  as  to 
the  number  of  days  required  to  dry  to  a  moisture  content  for 
safe  storage.  Fan  energy  requirements  (the  only  energy 
source  in  case  of  unheated  air)  are  significant  because  the 
fan  has  to  operate  for  much  longer  periods  of  time  than  with 
low  or  high  temperature  drying. 


2.4  COMBINATION  (HIGH-LOW  TEMPERATURE)  DRYING 

Combination  drying  is  a  system  in  which  high 
temperature  drying  is  followed  by  low  temperature,  or  on 
occasion,  ambient  air  drying.  The  purpose  of  high 
temperature  drying  is  to  reduce  the  grain  moisture  content 
to  a  level  where  drying  can  be  completed  without  the  risk  of 
spoilage  by  low  temperature  and  ambient  methods.  This  system 
of  drying  has  been  used  successfully  in  the  corn  belt  in  the 
United  States.  The  advantages  of  combination  drying  are 
reduced  energy  requirements,  increased  drying  capacity  of 
the  high  temperature  drier  and  improved  grain  quality  (Morey 
et  al.  1977;  Morey  et  al.  1978B). 

Morey  et  al.  (1978C)  conducted  an  extensive  field  scale 
evaluation  of  the  combination  drying  system  and  found  that 
the  energy  requirements  range  from  2100  to  3300  kJ/kg.  In  a 
similar  test,  Morey  et  al.  (1980)  reported  energy 
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requirements  in  the  range  of  1800  to  3800  KJ/kg;  the  higher 
statistic  occured  when  more  moisture  was  removed  in  the  high 
temperature  phase  than  in  the  low  temperature  phase.  The 
energy  requirements  are  considerably  less  than  the  hot  air 
driers,  moderately  higher  than  the  ambient  air  drying  and 
about  the  same  for  low  temperature  drying.  As  noted 
previously,  the  heat  energy  for  vaporization  comes  largely 
from  the  sensible  heat  of  the  ambient  air  in  the  low 
temperature  drying  phase  which  reduces  the  total  energy 
requirements  compared  to  high  temperature  drying. 

As  was  noted  above  the  combination  drying  has  been  used 
successfully  to  dry  corn  but  it  may  not  be  as  advantageous 
for  wheat  because  normally  the  latter  has  a  lower  moisture 
content  when  harvested. 


2.5  SOLAR  GRAIN  DRYING 

Grain  drying  is  a  process  which  can  be  easily  adapted 
for  the  use  of  solar  energy  (Bauman  and  Finner  1979).  With 
the  uncertainty  of  future  supplies  and  the  increases  in  the 
cost  of  fossil  fuels,  solar  energy  is  of  interest  as  a 
source  of  energy  for  grain  drying  (Feddes  et  al.  1980).  One 
major  factor  considered  in  the  installation  of  a  solar 
drying  system  is  the  capital  cost.  For  example,  Johnson  and 
Otten  (1980)  found  that  the  fixed  costs  of  the  collector 
were  sufficiently  high  to  make  solar  assisted  drying  more 
expensive  than  both  ambient  drying  and  conventional  high 
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temperature  drying  even  though  the  energy  cost  was  lowest 
for  solar  drying.  McLean  (1980)  reported  that  while  it  is 
extremely  pertinent  that  research  should  be  carried  out  in 
anticipation  of  rising  fuel  costs,  at  the  present  time  the 
installation  of  solar  panels  for  use  in  conjunction  with 
storage  drying  systems  is  uneconomical.  Hall  and  Davis 
(1979)  expressed  similar  views  and  stated  that  solar  energy 
systems  have  the  potential  for  minimizing  energy 
requirements  but  their  economic  justification  at  the  time  is 
margi na 1  . 

Peart  et  al.  (1980)  reported  fan  energy  requirements  of 
1800  to  2400  KJ/kg  for  solar  drying  of  corn  in  Indiana. 
Though  these  energy  requirements  are  considerably  less  than 
the  commonly  used  hot  air  driers  they  are  approximately  the 
same  as  ambient  air  drying.  The  reason  may  be  due  to  the 
solar  energy  gain  being  offset  by  the  higher  friction  losses 
in  the  solar  panels  and  the  additional  ducts.  One  advantage 
of  solar  drying  over  ambient  air  drying  is  a  higher  drying 
rate  and,  therefore,  the  potential  of  solar  drying  should 
not  be  neglected  in  future  research  on  grain  drying. 


2.6  MICROWAVE  VACUUM  DRYING 

Vacuum  drying  has  been  used  for  dehydrating  some 
agricultural  products  but  is  usually  restricted  to  foods 
where  quality  is  important  and;  therefore,  must  be  dried  at 
relatively  low  temperatures .  As  for  drying  cereal  grains 


under  vacuum  there  is  no  data  available  with  the  exception 
of  MIVAC,  the  experimental  microwave  vacuum  drying  system 
developed  recently  by  McDonnell  Douglas  Corporation  in  the 
United  States  (Gardner  and  Butler  1980). 

In  MIVAC  drying,  high  vacuum  and  microwave  energy  are 
used  to  dry  materials  quickly  and  at  temperatures  somewhat 
less  than  the  high  temperature  driers.  The  material  to  be 
dried  is  loaded  into  a  drying  chamber  which  is  then  isolated 
from  the  ambient  pressure  environment.  The  chamber  is 
evacuated  to  the  desired  pressure,  most  commonly  3.4  to  6.6 
KPa  (standard  atmospheric  pressure  is  about  101.5  kPa). 
During  drying,  the  material  flows  from  an  upper  hopper 
through  the  drying  zone  under  vacuum,  and  into  the  lower 
hopper  where  the  dried  material  is  collected. 

A  number  of  experiments  were  conducted  by  Gardner  and 
Butler  (1980)  over  a  sixteen  month  period  to  dry  many 
different  crops  such  as  corn,  peanuts,  sorghum,  rye,  and 
rice  using  MIVAC  dryer.  The  energy  requirements  were  in  the 
order  of  600  to  1000  KJ/kg  of  moisture  removed  which  are 
considerably  less  than  the  conventional  hot  air  driers  and 
somewhat  less  than  any  other  drying  system.  Another 
microwave  grain  drier  has  been  introduced  by  Ken  Bratney 
Company  of  Des  Moines,  Iowa  (Stauffer).  This  drier  basically 
works  on  the  same  principles  as  the  MIVAC  system,  and 
according  to  the  manufaturer,  it  uses  about  125  to  475  kJ/kg 
of  moisture  removed  from  soybeans.  The  commercial  success  of 
the  drier  is  assured  if  the  energy  requirements  are  as 
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c 1  a i med . 

The  energy  requirements  for  crop  drying  using  microwave 
vacuum  methods  appears  considerably  less  compared  to  other 
drying  methods  (as  quoted  by  the  two  manufacturers); 
however,  one  major  factor  that  a  farmer  considers  is  the 
initial  cost.  As  reported  by  BaKKer - Arkema  et  al.  (1978), 
the  initial  cost  is  not  competitive  with  conventional 
driers.  Boulanger  et  al.  (1969)  suggests  that  the  cost  of  a 
farm  sized  unit  which  can  process  7  cubic  meters  of  grain 
per  hour  (200  bu/hr)  is  in  the  order  of  $150,000  to 
$350,000.  Another  important  factor  limiting  the  use  of 
microwave  vacuum  drying,  is  the  use  of  electric  energy  for 
its  megnetron  unit.  It  has  been  shown  by  Morrison  (1978) 
that  electricity  is  still  a  relatively  expensive  form  of 
energy;  at  least  twice  the  equivalent  of  heat  energy  of 
fossil  fuels.  Also  a  skilled  operator  and  good  management 
appear  to  be  required  to  operate  this  system.  Due  to  these 
factors  and  the  possibility  of  high  maintenance  costs,  the 
use  of  these  driers  on  individual  farms  seems  unlikely  in 
the  near  future  unless  the  efficiency  is  extremely  good. 
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2.7  SUMMARY  AND  OBJECTIVE 

The  energy  requirements  for  a  number  of  drying  systems 
have  been  reported  and  discussed.  They  are  : 


1 . 

High  Temperature 

5000  (kJ/kg  water 

2. 

Low  Temperature 

3700-4400 

3. 

Combi nat ion 

1800-3800 

4. 

Sol  ar 

1800-2400 

5. 

Ambient 

900-2400 

6. 

Microwave  Vacuum 

600-1000 

The  energy  requirements  for  high  temperature  grain 
drying  are  almost  double  the  latent  heat  of  vaporization  of 
water  and  are  higher  than  any  other  system  in  use  whereas 
the  energy  requirements  for  vacuum  drying  system  using 
microwaves  were  the  lowest.  Largely  because  of  the  capital 
costs,  the  focus  will  likely  be  on  improving  the  efficiency 
of  the  other  systems. 

The  previous  discussions  show  that  the  energy 
requirements  for  ambient  air  drying  are  about  same  or  less 
than  the  latent  heat  of  water  vaporization.  As  mentioned 
earlier,  this  method  is  generally  considered  efficient  as 
the  heat  for  water  vaporization  comes  largely  from  the 
ambient  atmosphere.  One  major  problem  with  this  method  is 
the  lower  drying  rate  which  causes  a  higher  risk  of 
spoi 1  age . 

The  ambient  air  systems  usually  have  deep  beds  of  grain 
which  results  in  a  significant  pressure  gradient  across  the 
grain.  A  fan  is  used  to  provide  the  pressure  so  that  air 
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enters  the  bin  at  pressures  higher  than  atmospheric  and 
leaves  at  atmospheric  pressure.  Drying  occurs  on  the  basis 
of  the  ambient  air  conditions  at  atmospheric  pressure.  Since 
decreasing  the  total  pressure  will  increase  the  vapor 
pressure  differential  and  since  the  air  movement  through 
grain  is  based  only  on  the  static  pressure  difference  across 
the  grain,  one  would  expect  that  the  drying  capacity  of  an 
ambient  air  system  could  be  increased  without  affecting  the 
energy  input  by  placing  the  fan  on  the  outlet  of  the  grain 
bin  rather  than  at  the  inlet.  In  such  a  system,  the  pressure 
inside  the  grain  bin  will  be  slightly  below  atmospheric  (for 
example  in  the  order  of  85  to  93  KPa).  The  heat  energy  for 
water  vaporization  will  be  provided  by  the  ambient  air.  If 
the  energy  requirements  are  comparable  to  the  traditional 
ambient  drying  system  then  this  system  may  prove  a  better 
alternative;  consequently,  the  potential  of  this  system 
should  be  investigated. 

Based  on  the  above  discussions,  the  objective  of  this 
project  was  to  study  the  feasibility  of  drying  grain  with 
air  at  ambient  temperature  but  with  pressure  less  than 
atmospheric  and  compare  this  to  the  traditional  ambient  air 
drying  method.  The  basis  for  comparison  is  the  energy 
required  to  evaporate  one  Kilogram  of  water. 


3.  PROCEDURE 


3.1  INTRODUCTION 

Two  ambient  drying  systems  (above  and  below  atmospheric 
pressure)  were  studied  for  the  common  drying  period  in 
Western  Canada  of  August  through  November.  The  grain 
considered  was  wheat  because  of  the  large  quantity  grown  in 
Western  Canada.  The  study  method  was  to  calculate  the  energy 
requirements  for  removing  moisture  from  the  grain  by  the  two 
systems.  The  dry-bulb  and  wet-bulb  temperatures  (t  and  tw 
respectively)  prevailing  near  Edmonton,  Alberta,  as  given  by 
Cudbird  (1964),  were  used  for  the  study  (Appendix  A).  These 
temperatures  were  monthly  averages  for  every  hour  of  each 
day  for  the  four  months  mentioned  above. 

Seven  different  pressures  were  considered  for  this 
study,  with  four  below  and  three  above  atmospheric.  These 
pressures  were  86.7,  88.4,  90.1,  91.8,  95.2,  96.9,  and  98.5 
KPa  while  the  atmospheric  pressure  as  given  by  Cudbird 
(1964)  was  93.5  KPa.  Three  below  atmospheric  pressures 
correspond  to  three  of  the  four  above  atmospheric  pressures; 
having  same  magnitude  above  and  below  the  atmospheric 
pressure.  They  were  selected  to  compare  the  two  systems. 
These  pressures  represent  a  pressure  drop  of  1.7,  3.4,  5.1, 
and  6.8  KPa  through  grain  having  depths  of  approximately 
1.0,  1.3,  1.6,  and  1.9  meters  respectively  when  an  airflow 
of  3  m3/min- tonne  is  used  ( HuK ill  and  Shedd  1955).  A 
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pressure  drop  of  6.8  kPa  above  atmospheric  was  not 
considered  practical.  The  energy  requirements  for  a  system 
operating  at  atmospheric  pressure  are  not  included  because 
any  practical  drying  system  must  have  some  pressure  loss  in 
order  to  have  air  movement. 

Eight  initial  moisture  contents  were  selected  for  each 
moisture  content  from  15%  to  22%  and  the  final  moisture 
content  of  14%.  These  values  cover  the  most  likely  range  of 
moisture  contents  at  harvest  in  Western  Canada. 


3.2  ASSUMPTIONS 

A  number  of  assumptions  were  necessary  to  calculate  the 

energy  requirements.  They  are  : 

1)  Expansion  for  the  below  atmospheric  or 
exhausting  system  and  compression  for 
the  above  atmospheric  or  pressurizing 
system  is  isothermal;  that  is,  sufficient 
heat  is  gained  or  lost  to  the  surroundings 
to  maintain  the  temperature  in  the  plenum 
at  the  ambient  temperature  (see  figure  1). 

2)  No  change  in  enthalpy  occurs  as  the  air 
passes  through  grain. 

3)  The  airflow  is  low  enough  that  the  air 
leaves  grain  at  the  equilibrium  relative 
humidi ty  ( ERH ) . 
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Fan  location  ^Air  Out  ^ 


Figure  1:  Schematic  diagram  representing  the  two  drying 


systems . 
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4)  There  is  a  pressure  gradiant  across  the 

grain  mass . 

5)  The  fan  is  adiabatic. 

3.3  ENERGY  CALCULATIONS 

The  energy  requirements  were  calculated  on  the  basis  of 
water  removed  from  grain.  The  fan  energy  (E)  required  to 
move  air  through  grain  was  defined  as  kilojoule  per  kilogram 
of  dry  air  (kJ/kg  dry  air).  The  difference  in  the  humidity 
ratio  (kg  water/kg  dry  air)  of  the  air  entering  and  leaving 
the  grain  ( AX )  is  the  amount  of  moisture  removed  by  the  air 
from  the  grain  (figure  2)  as  noted  by  Harrison  (1969).  The 
figure  indicates  the  relationship  between  AX  and  the  air 
temperature  entering  and  leaving  the  grain  which  depends  on 
the  ERH.  Dividing  the  two  quantities  (E/aX)  would  yield  the 
energy  requirements  in  terms  of  heat  energy  per  unit  of 
moisture  removed  (kJ/kg  water  removed). 

A  computer  program  was  developed  to  calculate  the 
various  psychrometr i c  properties  of  air  before  and  after 
drying.  Calculations  for  the  humidity  ratio  of  the  air 
leaving  grain  (X7  )  required  the  use  of  an  iterative 
procedure.  A  complete  listing  of  the  program  is  provided  in 
the  Appendices  D  and  E. 

The  change  in  the  humidity  ratio  of  the  air  during 
drying  (aX)  was  calculated  for  various  moisture  content 
ranges  (22-14%,  21-14,  etc.  through  15-14%)  for  the  four 
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X' 


x 


Dry-Bulb  Temperature  (°C) 


Figure  2:  Psychrometr i c  Chart 


Humidity  Ratio  (kg  water/kg  dry  air) 
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drying  months  by  using  a  procedure  suggested  by  Harrison 
(1969).  The  dry-bulb  and  wet-bulb  temperatures  (t  and  tw 
respectively)  used  for  these  calculations  were  monthly 
averages  as  given  by  Cudbird  (1964).  If  aX  was  positive,  the 
air  had  a  drying  potential;  otherwise  it  had  a  wetting 
potential.  The  moisture  removal  per  day  was  determined  by 
dividing  the  sum  of  positive  humidity  ratios  by  the  number 
of  hours  the  changes  were  positive,  and  repeating  this  for 
moisture  content  changes  of  22  to  21%,  21  to  20%,  and  so  on 
until  15  to  14%.  For  the  moisture  range  16  to  14%  the 
moisture  removal  ( aX )  for  16  to  15%  was  added  to  the 
moisture  removal  for  15  to  14%.  For  17  to  14%  the  moisture 
removal  for  17  to  16%  was  added  to  the  moisture  removal  for 
16  to  14%  and  repeating  in  a  similar  manner  to  obtain  the 
moisture  removal  for  all  the  other  ranges. 

The  change  in  energy  of  the  air  as  it  passes  through  an 
adiabatic  fan  was  calculated  for  the  four  pressures  below 
atmospheric  and  three  above  as  mentioned  earlier.  The  air 
temperature  used  for  these  calculations  was  the  average  dry 
bulb  temperature  of  all  the  hours  when  aX  is  positive. 

The  procedure  for  calculating  aX  and  E  is  noted  below 
and  is  presented  in  a  flow  chart  in  Appendix  G  : 

STEP  1 

The  equilibrium  relative  humidity  ERH  (%)  for 
wheat  was  calculated  from  Chung's  equation 
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(Pfost  et  al.  1976;  Appendix  B)  : 

ERH  =  exp ( ( -  A/ ( t+C ) )  exp(-B  M))  (2) 

Where  M  =  grain  moisture  content  (%), 
t  =  dry-bulb  temperature  (°C), 

A  =  529.43  (hard  wheat), 

B  =  17.609,  and 
C  =  50.998. 

STEP  2 

The  ERH  values  were  determined  using  the  dry- 
bulb  temperatures  prevailing  near  Edmonton  as 
given  by  Cudbird  (1964)  and  the  hard  wheat  at 
nine  moisture  contents  of  14  through  22%. 

The  array  of  ERH  values  was  4  x  24  x  9. 

The  ERH  values  along  with  the  dry  and 
wet-bulb  temperatures  were  stored  in  a 
computer  file  (Appendix  C). 

STEP  3 

The  saturation  vapor  pressure  of  the  air 
Pws  (atmospheres)  over  ice  for  the 
temperature  range  of  -100°C  to  0°C 
was  calculated  from  (ASHRAE  1977)  : 


log(Pws)  =  -9 . 096936 ( 0  -  1 )- 3 . 56654  log ( 0 ) 

+  0 . 8768 1 7 ( 1  - ( 1/0) )-2. 2195983) 


(3) 
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and  for  the  temperature  range  of  0°C  to 
100°C  was  calculated  from  (ASHRAE  1977)  : 

log ( Pws )  =  1 0 . 79586 ( 1  -  0 ) +5 . 02808  log (  0) 

+  0.00015(1  -  10"8'29692((1/e)_1)) 

+  0 . 00042873  (  1 04-76955(1'e)  -  i  ) 

-2.219598  (4) 


Where  0  =  273 . 1 6/ ( tw+273 . 1 6 ) ,  and 

tw  =  wet-bulb  temperature  (°C), 
( Cudbi rd  1 964 ) . 


STEP  4 

The  humidity  ratio  of  the  air  at  saturation 
Xs  (lb  water/lb  dry  air)  was  calculated 
from  (ASHRAE  1977)  : 

Xs  =  0 . 62 1 98Pws/ ( B  -  Pws)  (5) 

and  B  =  (Bf  +  Batm)/2 

Where  B  =  mean  pressure  in  grain  mass  (KPa), 

Bf  =  one  of  the  seven  fan 
pressures  (KPa),  and 
Batm  =  atmospheric  pressure  (93.5  KPa). 
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STEP  5 

The  enthalpy  of  the  air  H  (kd/kg)  was 
calculated  from  (McQuiston  and  Parker  1977)  : 

H  =  tw  +  Xs ( 2503  .  1  +  1 . 86tw)  (6) 

STEP  6 

The  humidity  ratio  of  the  air  X  (lb  water/lb 
dry  air)  was  calculated  by  rearranging 
equation  (5)  and  also  replacing  tw  by  t. 

X  =  (H  -  t ) / ( 2503 . 1  +  1 . 86 t )  (7) 

The  humidity  ratio  before  and  after  compression 
or  expansion  does  not  change. 

STEP  7 

The  heat  gained  by  the  grain  or  change  in 
enthalpy  of  the  air  during  the  isothermal 
expansion  Q  (kd/kg),  was  calculated 
from  (Moyer  et  al .  1941)  : 

Q  =  R  ( t  +  273  .  16)  7  r?  (  P 1  /  P  2  )  (8) 

Where  P 1  , P 2  =  pressures  at  the  inlet  and 

exhaust  ( kPa ) ,  and 
R  =  universal  gas  constant. 
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Therefore,  the  enthalpy  of  air  (kJ/kg  dry  air) 
after  the  isothermal  expansion  is  (H'  )  : 

H'  =  H  +  Q  (9) 

There  was  no  change  in  the  enthalpy  of  the 
air  for  compression  because  the  heat  was 
assumed  to  be  lost  to  the  sur roundi ngs ; 
that  is,  isothermal  compression. 

STEP  8 

The  humidity  ratio  X'  (kg  water/kg  dry  air) 
of  the  air  leaving  grain  was  calculated  as 
fol lows  : 


i  ) 

calculate 

tout  = • H' /Ca 

where  Ca 

=  speci f i c  heat  of 

dry  air  (kj/kg  °  C ) , 

tout 

=  temperature  of  air 

leaving  grain  (‘C) . 

ii  ) 

decrement 

tout  by  TC 

iii  ) 

calculate 

Pws  at  tout  per 

equat i on 

(  3 )  or  ( 4  ) 

iv) 

calculate 

Xs  per  equation  (5) 

v) 

calculate 

H  per  equation  (6) 

Vi  ) 

calculate 

X'  per  equation  (7) 

vi  i  ) 

calculate 

the  relative  humidity 
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RH  (%)  by  (ASHRAE  1977)  : 

RH  =  1 00 ( X' /Xs ) / ( 1  - ( 1  - ( X'  /Xs ) Pws/B ) )  (10) 

viii)  go  to  step  (ii)  and  repeat  steps  (ii) 
to  (vii)  until  RH  is  equal  to  ERH. 

The  dry-bulb  temperature  of  step  (ii) 
and  the  humidity  ratio  of  step  (vi )  in 
the  final  loop  is  the  temperature  and 
the  humidity  ratio  (X'  )  of  the  air 
leaving  grain. 


STEP  9 

Calculate  aX  or  ( X ' - X )  then  go  to  step  (ii) 
and  repeat  steps  (ii)  to  (viii)  for  all 
the  nine  ERH  values  at  this  dry-bulb 
temperature.  Each  ERH  value  corresponds  to  a 
equilibrium  moisture  content  (Appendix  C). 

STEP  10 

Repeat  steps  (3)  to  (9)  for  all  the  dry-bulb 
temperatures  to  fill  the  array  (4  x  24  x  9). 

STEP  11 

Repeat  steps  (3)  to  (10)  for  all  the 
seven  pressures  to  fill  the  array 
(4  x  24  x  9  x  7) . 
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STEP  12 

Sum  the  positive  AX  values  for  each  hourly 
temperature  for  the  various  initial  moisture 
contents  of  22%,  21%,  until  15%.  The  hours 
(T)  when  AX  is  positive  were  also  summed. 
Dividing  s AX  by  T,  the  average  value  of 
AX  for  the  day  when  it  is  positive  is 
obtained . 

STEP  13 

Calculate  the  AX  values  for  different  initial 
moisture  content  ranges  of  22  to  14%, 

21  to  14%  through  15  to  14%  by  accumulating 
the  values  of  the  drying  potential  as 
noted  previously. 

STEP  14 

Determine  the  number  of  hours  (T)  available  to 
dry  each  day  for  different  moisture  content 
ranges  of  22-14%,  21-14%,  through  15-14%  by 
using  the  cummulative  values  of  T.  The 
procedure  followed  was  similar  to  the  one 
for  aX . 

STEP  15 

The  number  of  days  (D)  required  to  dry  grain 
from  various  initial  moisture  contents  to 
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14%  were  calculated  by  (Feddes  et  al.  1980; 
Harrison  1969)  : 


D  =  V  Ww  AW/q/T  AX 


(11) 


and  W  =  Mw  -  Mw'  ( 1 -Mw) / ( 1 -Mw'  ) 


(12) 


Where  V  =  specific  volume  of  dry  air 
(0.81  m3/kg ) , 

Ww  =  weight  of  grain  (1000  kg/tonne), 

o 

q  =  airflow  rate  (3  m  /mi n- tonne ) , 

T aX  =  kg  moi sture-mi n/kg  dry  air-day, 

AW  =  amount  of  moisture  removed  (kg), 
Mw  =  initial  moisture  content  (%),  and 
Mw' =  final  moisture  content  (%). 

These  equations  are  used  in  the  computer 
program  shown  in  Appendix  F. 


STEP  16 

The  change  in  energy  of  the  air  E  (kJ/kg 
dry  air)  after  it  passes  through  an  adiabatic 
fan  was  calculated  from  ( De  Nevers  1971)  : 


E  =  ( R ( t  +  273 . 1 6 ) k/ ( k- 1 ) ) ( (P2/P1 ) (k-1)  /k - 1  )  (13) 


Where  P 1 , P 2  =  pressures  at  the  inlet  and 

exhaust  ( kPa ) , 
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R  =  universal  gas  constant,  and 
K  =  1  .4 


STEP  17 

The  energy  requirements  for  drying  the  grain 
(kd/kg  water  removed)  were  calculated  by 
dividing  E  by  aX . 


4.  DISCUSSION  OF  RESULTS 


The  energy  required  to  remove  one  hi logram  of  moisture, 
the  number  of  hours  per  day  when  drying  can  occur,  and  the 
number  of  days  required  to  dry  wheat  to  14%  are  shown  in 
tables  1,  2  and  3.  Table  1  shows  that  some  of  the  energy 
requirements  per  Kilogram  of  moisture  are  less  than  the 
latent  heat  of  vaporization  for  grain  moisture  (2800  KJ/Kg) 
as  given  by  Brooker  et  al.  (1974).  As  noted  in  chapter  2, 
much  of  the  latent  heat  for  vaporization  comes  from  ambient 
air  and  is  not  included  in  the  energy  calculations.  Table  2 
shows  that  more  hours  per  day  are  available  for  drying  when 
the  pressure  is  below  atmospheric  than  above.  In  addition, 
drying  to  14%  on  a  24  hour  basis  is  possible  for  any  of  the 
four  drying  months  if  the  pressure  is  88.4  KPa  or  less.  This 
is  particularly  significant  for  the  month  of  November  when 
grain  cannot  be  dried  to  14%  in  the  conventional  manner 
(above  atmospheric  pressure)  using  ambient  air.  Table  3 
indicates  that  the  number  of  days  required  to  dry  wheat  from 
various  initial  moisture  contents  to  14%  are  generally  less 
for  pressures  below  atmospheric  than  above.  Also  table  3 
indicates  that  drying  to  14%  cannot  be  completed  during 
October  at  22%  initial  moisture  content  and  at  any  initial 
moisture  content  during  November  for  the  above  atmospheric 
pressures . 

The  energy  requirements  for  the  months  of  August, 
September,  and  October  at  three  initial  moisture  contents  of 
15,  18,  and  21%  are  shown  graphically  in  figures  3,  4,  and 
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Drying  to  14%  moisture  content  cannot  be  achieved  at  the  above  atmospheric  pressures. 
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Drying  to  14%  moisture  content  cannot  be  achieved  at  the  above  atmospheric  pressures. 
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*  Drying  to  14%  moisture  content  cannot  be  achieved  at  the  above  atmospheric  pressures . 
**  Drying  to  14%  moisture  content  cannot  be  achieved  during  one  month  period. 
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Figure  3:  Effect  of  Pressure  and  Initial  Moisture  Content 
on  the  Energy  Requirements  during  August 
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Figure  4:  Effect  of  Pressure  and  Initial  Moisture  Content 
on  the  Energy  Requirements  during  September 
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Figure  5:  Effect  of  Pressure  and  Initial  Moisture  Content 
on  the  Energy  Requirements  during  October 
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5.  November  is  not  included  since  drying  to  14%  cannot  be 
achieved  during  November  at  pressures  above  atmospheric. 
These  figures  show  that  the  energy  requirements  increase  as 
the  pressure  drop  through  the  grain  increases  for  both  the 
ambient  air  systems  evaluated.  The  increases  for  either 
pressurizing  or  exhausting  (above  and  below  atmospheric)  are 
due  to  the  effort  to  pump  the  air  through  grain;  that  is, 
the  energy  requirements  for  the  ambient  air  system  are 
attributed  to  the  friction  loss  of  the  air  through  the 
grain.  Increasing  the  grain  depth  or  the  airflow  or  both 
would  result  in  increased  energy  requirements.  The  rate  of 
increase;  however,  is  greater  for  pressures  above 
atmospheric  than  below.  In  other  words,  the  greater  the 
depth  of  grain  or  airflow,  the  higher  the  energy 
requirements  for  either  exhausting  or  pressurizing  the  bin 
but  the  rate  of  increase  is  less  for  exhausting  than  for 
pressurizing.  The  reason  for  this  is  noted  later. 

The  energy  requirements  increase  as  the  initial 
moisture  content  of  grain  decreases1  as  shown  in  figures  3, 
4,  and  5.  This  relationship  is  attributed  to  the  variation 
in  the  bonding  forces  within  the  Kernel  holding  the 
moisture.  For  wet  grain  the  water  vaporizes  as  it  does  from 
a  liquid  surface  but  for  dry  grain  the  water  is  held  with 
considerable  force.  This  increase  was  observed  for  both  the 
above  and  below  atmospheric  pressures;  however,  the  increase 

1  The  energy  to  dry  given  quantity  of  grain  increases  with  an 
increase  in  initial  moisture  content. 
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for  pressures  above  atmospheric  was  higher  than  for 
pressures  below  atmospheric. 

The  energy  requirements  increase  as  the  drying  season 
proceeds,  mainly  due  to  the  drop  in  ambient  air  temperature. 
Average  daily  air  temperatures  drop  from  about  15.5  to  4.4°C 
at  Edmonton  during  the  period  August  to  October  thereby 
reducing  the  drying  potential  of  the  air.  The  increase  in 
energy  requirements  is  slightly  greater  for  pressures  above 
atmospheric  than  below  (figures  3,  4,  and  5).  For  example,  a 
comparison  of  figures  3  and  5  indicates  that  the  energy 
requirements  increased  50%  for  the  period  August  to  October 
for  a  pressure  above  atmospheric  against  a  40%  increase  for 
the  same  pressure  below  atmospheric. 

More  important  than  the  above  is  the  indication  from 
figures  3,  4,  and  5  that  less  energy  per  Kilogram  of 
moisture  removed  is  required  by  operating  at  pressures  below 
atmospheric  than  above.  This  can  be  further  illustrated  in 
figures  6,  7,  and  8  where  efficiency  gain  was  used  to 
express  the  advantage  in  using  pressures  less  than 
atmospheric.  In  these  figures  the  pressures  above  and  below 
atmospheric  are  replaced  with  appropriate  depths  of  grain  as 
indicated  in  chapter  3  namely  1.0,  1.3,  and  1.6  meters.  The 
efficiency  gain  is  defined  as  the  difference  in  energy 
requirements  between  the  two  systems  expressed  as  a 
percentage  of  the  energy  requirements  of  the  system  using 
pressures  above  atmospheric.  The  figures  show  that 
efficiency  gain  increases  with  the  depth  of  grain.  This 
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Gain  in  Drying  Efficiency  by  exhausting  the  bin 
rather  than  pressurizing  during  August 
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Gain  in  Drying  Efficiency  by  exhausting  the  bin 
rather  than  pressurizing  during  September 
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suggests  that  using  the  fan  to  exhaust  a  bin  type  drier  is 
more  efficient  than  using  the  fan  to  pressurize  it,  and  the 
efficiency  gain  increases  with  the  depth  of  grain.  The  gain 
is  attributed  to  the  following  : 

1.  Heat  is  acquired  from  the  surroundings  during  the 
isothermal  expansion  (below  atmospheric)  and  is  lost  to 
the  surroundings  during  isothermal  compression  (above 
atmospheric);  therefore,  the  heat  energy  for  exhausting 
the  bin  is  less  than  for  pressurizing  it. 

2.  A  greater  vapor  pressure  differential  is  created 
between  the  grain  and  its  surroundings  when  the  total 
pressure  is  reduced  below  the  atmospheric  pressure  and 
due  to  this  greater  vapor  pressure  differential,  the 
difference  in  the  humidity  ratios  of  the  air  entering 
and  the  air  leaving  grain  is  greater;  therefore,  more 
moisture  can  be  removed  from  the  grain  with  the  same 
amount  of  air. 

3.  The  latent  heat  of  vaporization  is  lower  for 
pressures  below  atmospheric  than  above;  therefore, 
requiring  less  energy. 


As  stated  previously,  the  energy  requirements  increase 
at  a  greater  rate  when  pressures  are  increased  above 
atmospheric  than  when  decreased  an  equal  amount  below 
atmospheric.  This  may  be  explained  by  examining  each  system 
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separately.  In  case  of  the  above  atmospheric  pressure  system 
when  the  pressure  is  increased  (either  by  increasing  grain 
depth  or  the  airflow),  the  resistance  offered  by  the  grain 
increases  greatly;  which  causes  the  energy  requirements  to 
increase  at  a  substantial  rate.  With  the  below  atmospheric 
system  when  the  pressure  is  further  decreased,  the  energy 
requirements  to  operate  the  system  increase;  however,  the 
vapor  pressure  differential  created  also  increases.  This 
increase  in  vapor  pressure  differential  partially  offsets 
the  higher  energy  requirements  due  to  a  higher  airflow  or 
greater  grain  depth;  hence,  the  resulting  increase  in  the 
energy  requirements  for  exhausting  is  less  than  for 


pressur i zing . 


• 

5.  SUMMARY  AND  CONCLUSIONS 

The  energy  requirements  for  drying  wheat  during  the 
period  August  through  November  at  three  above  and  four  below 
atmospheric  pressures  for  the  ambient  air  conditions 
prevailing  near  Edmonton,  Alberta,  were  calculated.  The 
basis  for  calculating  the  energy  requirements  was  the  amount 
of  heat  energy  required  to  remove  one  Kilogram  of  water  from 
grain.  An  empirical  equation  from  Pfost  et  at  (1976)  and  a 
set  of  equations  from  ASHRAE  (1977)  were  used  respectively 
to  calculate  the  equilibrium  relative  humidity  and  the 
various  psychrometr i c  properties  of  the  air  entering  and 
leaving  grain.  The  energy  required  to  pump  the  air  through 
grain  was  calculated  by  assuming  an  airflow  and  several 
pressure  losses.  The  heat  energy  per  unit  of  moisture 
removed  (KJ/kg  water  removed)  was  determined  by  dividing  the 
energy  by  the  drying  capacity.  Based  on  these  calculations, 
the  following  conclusions  are  drawn  : 

1.  The  energy  requirements  for  drying  are  less  at 
pressures  below  atmospheric  than  above. 

2.  The  energy  requirements  increase  at  a  greater 

rate  for  pressures  above  atmospheric  than  below  as  the 
the  depth  of  grain  or  airflow  is  increased. 

3.  Drying  can  continue  (to  14%  moisture  content)  for 
more  hours  per  day,  and  even  on  a  24  hour  basis,  using 
pressures  below  atmospheric  for  the  usual  drying  period 
in  Western  Canada. 
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4.  Number  of  days  required  to  dry  wheat  to  14% 

moisture  content  are  generally  less  for  pressures  below 
atmospher i c . 

5.  The  energy  efficiency  is  greater  for  pressures  below 
atmospheric  than  above  and  this  differential  increases 
for  decreasing  initial  moisture  content. 

6.  The  energy  efficiency  decreases  as  the  drying  season 
proceeds  but  the  decline  proceeds  at  a  slightly  lower 
rate  for  pressures  below  atmospheric  than  above. 


6.  SUGGESTIONS  FOR  FURTHER  STUDY 


The  following  suggestions  could  prove  useful  in  further 
investigating  and  improving  the  drying  efficiency  of  the  low 
pressure  drying  system.  Firstly,  an  experimental  drier 
should  be  designed  and  tested  at  certain  above  and  below 
atmospheric  pressures  used  in  this  study;  such  as,  86.7, 
91.8,  95.2,  and  98.5  KPa.  Some  additional  pressures  may  be 
included  in  the  study  to  enhance  the  application  of  the 
results.  These  pressures  can  be  established  by  using 
different  grain  depths  and/or  airflow  rates.  The  grain  depth 
or  airflow  values  considered  should  be  the  one  commonly  used 
in  Western  Canada  for  the  conventional  manner  of  ambient  air 
drying;  that  is,  to  pressurize  the  grain  bin.  The 
exper imenta 1 -set  up  should  measure  the  energy  input  for 
moving  air  through  the  grain,  the  dry-bulb  and  dew-point 
temper atures ,  and  the  moisture  loss  from  the  grain.  The 
possibility  of  heating  the  air  a  few  degrees  above  ambient 
during  the  adverse  weather  conditions  which  increases  the 
drying  potentional  of  the  air  may  also  be  considered. 
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APPENDIX  A 


Hourly  Dry  and  Wet-Bulb  Temperatures  (F)  prevailing 
near  Edmonton,  Alberta  during  the  four  drying 
months  of  August  through  November,  as  given  by 
Cudbird  (1964)  for  24  hours  starting  at  midnight 

AUGUST 

D.B.  Temperatures  :  56.5,55.5,54.5,53.8,53.0,52.4, 

52.3,54.1,56.8,59.7,62.1,63.9,65.9,67.1,68.1,68.7, 

68.7,68.5,67.3,65.6,62.9,60.7,58.9,57.7 

W.B.  Temperatures  :  52.9,52.2,51.7,51.1,50.5,50.1, 

50.2,51.4,52.9,54.4,55.5,56.0,56.7,57.0,57.3,57.4, 

57.3,57.3,57.0,56.6,55.8,55.0,54.2,53.5 

SEPTEMBER 

D.B.  Temperatures  :  47.5,46.5,45.6,44.8,44.1,43.6, 

42.8.43.4.45.7.48.7.51.7.53.9.56.1.57.7.58.8.59.3, 
59.1,58.7,56.7,54.0,51.9,50.2,48.9,48.1 

W.B.  Temperatures  :  44.0,43.4,42.8,42.2,41.8,41.4, 

40.9.41.3.42.9.44.7.46.3.47.3.48.3.48.9.49.2.49.3, 
49.4,49.2,48.4,47.3,46.3,45.4,44.7,44.3 

OCTOBER 

D.B.  Temperatures  :  37.6,36.7,36.0,35.4,34.9,34.6, 

34.1,34.0,35.0,37.8,41.1,43.6,46.3,48.0,49.0,49.5, 

49.0,47.6,44.7,42.9,41.4,40.0,39.0,38.2 

W.B.  Temperatures  :  34.3,33.7,33.2,33.8,32.4,32.1, 

31.8,31.7,32.5,34.5,36.5,38.0,39.5,40.2,40.6,40.8, 

40.6,39.9,38.4,37.5,36.6,35.8,35.1,34.6 

NOVEMBER 

D.B.  Temperatures  :  23.9,23.5,23.1,22.9,22.6,22.4, 

22.0,22.0,22.0,22.9,24.9,26.7,29.2,30.7,31.5,31.4, 

30.4,28.9,27.4,26.5,25.6,25.0,24.5,24.1 

W.B.  Temperatures  :  22.0,21.7,21.4,21.2,20.9,20.8, 

20.6.20.4.20.4.21.1.22.5.23.9.25.7.26.7.27.3.27.3, 
26.6,25.7,24.7,24.0,23.4,22.9,22.4,22.1 
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APPENDIX  B 


C  The  program  ERH  was  used  to  calculate  the 
C  Equilibrium  Relative  Humidities  for  wheat 
C  under  the  ambient  temperatures  and  various 
C  initial  moisture  contents  as  noted  in  chapter  3 
C 

DIMENSION  A ( 24 ) , B ( 24 ) , T ( 24 ) 

REAL  MW ( 9  )  , MD ( 9 ) , ERH ( 9 ) 

INTEGER  I,J 
D  =  529 . 43 
B  =  1  7 . 609 
C  =  50 . 998 

RE AD ( 5  ,  10)  (A(I) ,1  =  1,24) , (B(I) ,1  =  1,24) 

10  FORMAT ( F5 . 3 ) 

READ  (  5  ,  1  1  )  (MW(vJ)  ,0=1,9) 

11  FORMAT ( 9F4 . 2 ) 

DO  30  1=1,24 
T(I)=(A(I)-32. ) *5 . /9 . 

DO  12  0=1,9 

MD ( 0 ) = 1 00 . *MW ( 0 ) / ( 100. -MW( 0 ) ) 

MD ( 0 ) =MD ( 0 ) / 1 00 . 

ERH(0)=EXP( (-D/(T(I)+C) ) *EXP ( -B*MD ( 0 ) )) *100.0 

12  CONTINUE 
TEMP= A ( I ) 

TEMP  1 =  B ( I ) 

WRITE (6, 20)  (TEMPI , TEMP , ( ERH ( 0 ) , 0= 1 , 9 ) ) 

20  FORMAT ( 1 1 F8 . 1 ) 

30  CONTINUE 
STOP 
END 
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APPENDIX  D 


C  The  program  PSYALL  performs  psychromet r i c 
C  calculations  during  the  drying  process  when 
C  pressure  surrounding  grain  is  below  atmospheric 
C 

LOGICAL* 1  LIS ( 1 )/' *' / 

DIMENSION  ERH ( 9  ) 

WR I TE ( 6 , 4 ) 

4  FORMAT ('  1  TIME  MOIS  BARO  ENTH  TWET  TIN  TOUT', 
I7  DELT  WIN  WOUT  DELW  ERH '/) 

WR I TE ( 7 , 96 ) 

96  FORMAT ( ' ENTER  THE  REDUCED  PRESSURE.') 

READ ( 7 , L I S  )  B2 
C  B2  =  ( Bf an  +  Batm)/2.0 
C  assume  linear  pressure  gradient 
DO  9999  J=1 ,24 
B  =  27 .6 

READ (5,3)  TW , T 1 , ( ERH ( I ) , 1=1 ,9) 

WR I TE ( 7 , 3  )  TW , T 1 , ( ERH ( I ) ,1  =  1 ,9) 

3  FORMAT ( F4 . 1 , 10F8. 1 ) 

COMPUTE 

CALL  PRESS ( T 1 ,PS) 

CALL  PRESS ( TW , PX ) 

WX=.622*PX/(B-PX) 

H= . 24*TW+WX* ( 1061 . +.45*TW) 

W 1  =  ( H- . 24*T 1 )/(  1061  .  +  (  .  45*T 1 ) ) 

Completed  the  first  step  -  At  this  stage,  we  have  inlet 
C  absolute  humidity  W1  and  enthalpy  H.  Now  decrease 
C  the  pressure.  The  enthalpy  of  air  will  change 

C  because  of  the  heat  gained  (Q)  from  the 

C  surroundings  due  to  isothermal  expansion. 

TH 1  =H/0 . 24 

Q=(53.3*(T 1+460. 0)*ALOG(B/B2) )/778. 

HBAR=H+Q 
TH2=HBAR/ 0 . 24 
DELH=TH2-TH 1 
TBAR=T 1 +DELH 

C  HBAR  and  TBAR  are  new  values  of  enthalpy 
C  and  temperature  of  air  at  low  pressure. 

C 

C  DO  LOOP  for  calculating  the  absolute  humidity 
C  of  air  (WOUT)  leaving  grain  at  the  ERH. 

DO  900  1=1,9 
TM=  HBAR/. 24 
X=1  .0 

21  TM=TM-X 

CALL  PRESS ( TM , PS ) 

W= ( HBAR - . 24*TM ) / ( 1061 .+( .45*TM) ) 

WS=.622*PS/ (B2-PS) 

C  CALL  RELHUM ( B2 , W , WS , PS , R 1 )  continued  .. 
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u=w/ws 

I F ( U . GT . 1 .  )  U=  1  . 

RH= 1 00 . *U/ ( 1 . - ( ( 1 . -U)*( PS/B2)  )  ) 

IF(RH.LT.ERH(I) )  GO  TO  21 
IF(X.EQ.O. 1 )  GO  TO  23 
TM=TM+X 
X=X/10 
GO  TO  21 
23  CONTINUE 

C  - OUTPUT - 

W=W*1000. 

W 1 = W 1  *  1 000 . 

DELW=W-W1 
DELT  =  TBAR - TM 
L  =  23- 1 

WRITE (6, 34)  J,L,B2,H,TW,T1 , TM, DELT ,W1 ,W,DFLW,ERH(I) 
34  FORMAT ( 1X,I2,5X,I2, IX,  10F10.2) 

W 1 = W 1 / 1 000 
900  CONTINUE 
9999  CONTINUE 
STOP 
END 

C  - Functions  and  subroutines - 

SUBROUTINE  PRESS ( T  ,  PP ) 

TT  =  273. 16/ ( ( ( T - 32 . ) *5 . /9  .  )+273. 16) 

I F ( T . LT . 32 .  )  GO  TO  1 

P= 1 0 . ** ( 10. 79586* ( 1 .-TT)+5.02808*ALOG10(TT) 

1  +  1 . 50474E-4* ( 1 . - 1 0** ( - 8 . 29692* ( (  1  . / T T  )  - 1  .  )  )  ) 

2  +  0 . 42873E -  3* ( 1 0 . ** ( 4 . 76955* ( 1 . -TT ) ) -  1 . ) -2.219598) 

GO  TO  2 

1  P= 1 0.**( -9. 096936* (TT-1 . ) - 3 . 56654*ALOG 1 0 ( TT  ) 

1+0 .8768 17* ( 1 . - ( 1 . /TT) ) -2.219598) 

2  PP  =  P*29 . 92 1 
RETURN 

END 


APPENDIX  E 


C  The  program  PSYALL1  performs  psychrometr i c 
C  calculations  during  the  drying  process  when 
C  pressure  surrounding  grain  is  above  atmospheric 
C 

LOG  I  CAL*  1  LIS(1)/'*'/ 

DIMENSION  ERH ( 9  ) 

WR I TE ( 6 , 4 ) 

4  FORMAT ( ' 1 T IME  MOIS  BARO  ENTH  TWET  TIN  TOUT', 
T  DELT  WIN  WOUT  DELW  ERH'/) 

WR I TE ( 7 , 96 ) 

96  FORMATS  ENTER  THE  INCREASED  PRESSURE.') 

READ ( 7 , L I S )  B2 
C  B2  =  ( B  f an  +  Batm)/2.0 
C  assume  linear  pressure  gradient 
DO  9999  J=  1  r  24 
B  =  27 . 6 

READ (5,3)  TW , T 1 , ( ERH ( I ) , 1=1 ,9) 

WR I TE ( 7 , 3  )  TW , T 1 , ( ERH ( I ) , 1  =  1 ,9) 

3  FORMAT ( F4 . 1 , 10F8. 1 ) 

COMPUTE 

CALL  PRESS ( T 1 , PS ) 

CALL  PRESS ( TW , PX ) 

WX= . 622*PX/ ( B- PX ) 

H= . 24*TW+WX* ( 1061 .+.45*TW) 

W 1  =  ( H- . 24*T 1 )/(  1061  .  +  ( . 45*T  1  )  ) 

Completed  the  first  step  -  At  this  stage,  we  have  inlet 
C  absolute  humidity  W1  and  enthalpy  H. 

C  Now  air  enters  grain  where  the  pressure  is 
C  above  atmospheric  and  W1  &  H  remains  unchanged. 

C 

C  DO  LOOP  for  calculating  the  absolute  humidity  WOUT 
C  and  temperature  TOUT  of  air  leaving  grain  at  ERH. 

C 

DO  900  1=1,9 
TM=  H/.24 
X=1  .0 

21  TM=TM-X 

CALL  PRESS ( TM , PS ) 

W=(H-.24*TM)/( 1061 .+( .45*TM) ) 

WS= . 622*PS/ ( B2- PS ) 

C  CALL  RELHUM(B2,W,WS,PS,R1 ) 

U=W/WS 

IF(U.GT. 1 . )  U= 1 . 

RH= 1 00 . *U/ ( 1 .-( ( 1 .-U)*(PS/B2) ) ) 

IF(RH.LT.ERH(I) )  GO  TO  21 
I F ( X . EQ . 0  .  1  )  GO  TO  23 
TM=TM+X 
X=X/10 

GO  TO  21  conti nued  . . 
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23  CONTINUE 

C  - OUTPUT - 

W=W* 1000 . 

W1=W1*1000. 

DELW=W- W 1 
DELT  =  T  1  -  TM 
L  =  23- 1 

WRITE (6, 34)  J,L,B2,H,TW,T 1 ,TM,DELT,W1 ,W,DELW,ERH( I ) 
34  FORMAT ( IX, I2,5X, 12, IX, 10F10.2) 

W 1 = W 1 / 1 000 
900  CONTINUE 
9999  CONTINUE 
STOP 
END 

C  - Functions  and  subroutines - 

SUBROUTINE  PRESS(T,PP) 

TT=273 . 16/ ( ( (T-32. )*5./9. )+273. 16) 

I F ( T . LT . 32 .  )  GO  TO  1 

P= 1 0 . ** ( 10. 79586* ( 1 .-TT)+5.02808*ALOG10(TT) 

1  +  1 . 50474E-4* ( 1 . - 10** ( -8. 29692* ( ( 1 . / T T  )  -  1  .  )  )  ) 

2  +  0 . 4287  3E -  3* ( 1 0 . ** ( 4 . 76955* ( 1 . -TT ) ) -  1 . ) -2.219598) 

GO  TO  2 

1  P=10.**( -9. 096936* (TT-  1 .  )  - 3 . 56654*ALOG 1 0 ( TT ) 

1+0. 87681 7* ( 1 .-( 1 ./TT) ) -2.219598) 

2  PP  =  P*29 . 92 1 
RETURN 

END 


APPENDIX  F 


C  The  program  FILE1  was  used  to  calculate  the  energy 
C  requirements  per  Kilogram  of  water  removed,  the 
C  number  of  hours  available  each  day  to  dry,  and  the 

C  number  of  days  required  to  dry  wheat  to  14%  as 

C  suggested  by  Feddes  et  al.  (1980). 

C 

DIMENSION  DELW( 8 , 24 ) , SUM  1(8) , SUM2 ( 8 ) ,SUM3(8)  , 

1  TEMP (8, 24) , SUM4 ( 9 ) , SUM5 ( 9 ) ,SUM6(9) ,DELWB(8) , 

1 TEMPB ( 8 ) ,DELWR(8) ,TEMPR(8) , TR ( 8 ) , P  2 ( 7 ) , AEN ( 8 ) , 

1 AKJ ( 8 ) , DA YS ( 8 ) , AA ( 8 ) , TDELWR ( 8 ) 

DATA  AA/0. 15,0. 16,0. 17,0. 18,0. 19,0.20,0.21 ,0.22/ 
DATA  P2/29 . 1,28.6,28.1,27.1,26.6,26.1,25.6/ 

DO  99  L=1  ,7 

C  I=INITIAL  MOISTURE  CONTENT 
C  K=H0UR 

DO  14  K=  1  ,24 
DO  13  1=1,8 

READ ( 5 , 1 1 )  TEMP ( I ,K) , DELW ( I , K ) 

C  WR I TE ( 7 , 97 )  TEMP ( I , K ) , DELW ( I , K ) 

11  FORMAT (43X, F 10.2,40X, F 10.2) 

C  97  FORMAT(43X,F10.2,40X,F10.2) 

13  CONTINUE 

14  CONTINUE 

DO  10  K=1 ,8 
SUM  1 ( K ) =0 . 0 
SUM2 ( K ) =0 . 0 
SUM3 ( K  )  =0 . 0 
SUM4 ( K ) =  0 . 0 
SUM5 ( K ) =0 . 0 
SUM6 ( K ) =0 . 0 
DO  20  1=1,24 

C  SUM1  is  sum  of  humidity  ratios  per  day  for  each 
C  K  initial  moisture  content 

C  SUM2  is  number  of  drying  hours  per  day  for  each 
C  K  initial  moisture  content 

C  SUM3  is  sum  of  dry-bulb  temperatures  when  the 
C  air  has  drying  potential 

IF (DELW(K, I )  . GT .  0.0)  SUM  1 ( K ) =  SUM  1 (K)+DELW(K, I ) 

I F ( DELW ( K , I ) . GT .  0.0)  SUM2 ( K ) =SUM2 ( K ) + 1 . 0 
I F ( DELW ( K , I ) . GT .  0.0)  SUM3 ( K ) =SUM3 ( K ) +TEMP ( K , I ) 

20  CONTINUE 

DE  LWB ( K ) =  SUM  1 (K)/SUM2(K) 

TEMPB(K)=SUM3(K)/SUM2(K) 

10  CONTINUE 

C 

DO  15  K=1 ,8 
KK=9-K 

SUM4 ( K+ 1 ) =SUM4 ( K ) +DELWB ( KK ) 

SUM5 ( K+ 1 ) =SUM5 (K)+TEMPB(KK)  conti nued 


60 


61 


SUM61K+1 ) =SUM6 ( K )  +SUM2 ( KK ) 

C  SUM4  is  the  accumulative  sum  of  the 
C  humidity  ratios(DELTA  X's) 

C  SUM5  is  the  accumulative  sum  of  the 
C  dry-bulb  temperatures 

C  SUM6  is  the  accumulative  sum  of  no  of  hours 
D I  V  =  K 
P  1  =27 . 6 
R=53 . 3 
A  K  = 1  .4 
Q=  3 . 0 

C  Q  the  airflow  rate  in  CFM  per  bushel 
DELWR ( K ) =  SUM4 ( K+ 1 )/DIV 
TEMPR ( K ) =SUM5 ( K+ 1 )/DIV 
TR  (  K  )  =  SUIV16  (  K+  1  )  / D I  V 
TDELWR(K)=TR(K)*DELWR(K)*60. 0/1000.0 
C  DELWR  is  the  average  DELTA  X  value  for  various  m.c  ranges 
C  TEMPR  is  the  average  D.B.TEMP  value  for  various  m.c  ranges 
C  TR  is  the  average  NO  of  HOURS  for  various  m.c  ranges 

C  TDELWR  is  the  product  of  DELWR  &  TR.  It's  units  are 
C  lb  water  x  min  /  lb  dry  air  per  day 

AEN ( K ) =ABS ( ( ( ( TEMPR (K)+460. ) *R*AK ) / ( AK- 1 . 0) )* 

1  ( (P2(L)/P1  )**( (AK-1 . 0 ) /AK ) -  1 .0) 1/778.0 
DELM  =  AA ( K ) -  0 . 1 4* ( 1 . 0-AA ( K ) ) /0 . 86 
DAYS ( K ) = ( 13.0*60.0*DELM)/ ( Q*TDELWR ( K ) ) 
AKJ(K)=(AEN(K)/DELWR(K) )*1000. *2. 326 
C  AEN  is  the  energy  required  to  move  air  through  grain 
C  AKJ  is  the  KJ/Kg  of  water  removed 
C  DAYS  is  the  number  of  days  to  dry  grain 
15  CONTINUE 

C  Now  output  all  values  for  this  month  and  pressure. 

WRITE (6, 12)  L , ( DELWR ( K 1 , K  = 1 , 8 1 , ( TEMPR (K)  ,K=1 ,8) , 

1  ( T  R ( K )  ,K=1 ,8) , (AEN(K) , K  = 1 ,8) , (AKJ(K) , K= 1 ,8) , 

1  ( DAYS ( K ) ,K=1 ,8) 

12  FORMAT (2X. I4,8F10.2,/8X,8F10.2,/8X,8F10.2, 

1/8X,8F 10.2, /8X,8F 10.2 , /8X,8F  10.2) 

99  CONTINUE 
STOP 
END. 


' 


APPENDIX  G 


Flow  chart  showing  in  sequence  the  different  computer 
programs  used  to  perform  various  calculations 
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